Entanglement concentration of continuous variable quantum states 
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We propose two probabilistic entanglement concentration schemes for a single copy of two-mode 
squeezed vacuum state. The first scheme is based on the off-resonant interaction of a Rydberg 
atom with the cavity field while the second setup involves the cross Kerr interaction, auxiliary mode 
prepared in a strong coherent state and a homodyne detection. We show that the continuous- 
variable entanglement concentration allows us to improve the fidelity of teleportation of coherent 
states. 
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I. INTRODUCTION 

Quantum entanglement is an essential ingredient of 
many protocols for quantum information processing such 
as quantum teleportation jl], |^] or quantum cryptography 
||. In order to achieve optimum performance of these 
protocols, the two involved parties, traditionally called 
Alice and Bob, should share a pure maximally entangled 
state. In practice, however, we are often able to generate 
only non-maximally entangled state and, in addition, the 
distribution of the entangled state between the two dis- 
tant parties via some noisy quantum channel will degrade 
the entanglement and Alice and Bob will share some par- 
tially entangled mixed state. One of the most important 
discoveries in the quantum information theory was the 
development of the entanglement purification protocols 
that allow Alice and Bob to extract a small number of 
highly entangled almost pure states from a large number 
of weakly entangled mixed states @, ||[ ^| . These proto- 
cols involve only local operations and classical communi- 
cation (LOCC) between the two parties, therefore they 
can be performed after the distribution of the entangled 
states. 

In the simplest scenario Alice and Bob share a 
pure non-maximally entangled state in a d-dimensional 
Hilbert space whose Schmidt decomposition reads 



c,- ai)A 



\0i 



(1) 



where each set of states \ctj) and \(3j) forms a basis. Al- 
ice and Bob would like to prepare from \yb) a state with 
higher entanglement by means of LOCC operations. Re- 
markably, this is possible, albeit only with certain prob- 
ability, even if they share only a single copy of this state. 
The procedure that accomplishes this task was fittingly 
called the Procrustean method § , because it cuts off the 
Schmidt coefficients Cj to the size of the smallest one. In 
this way, Alice and Bob obtain, with certain probability, 
a maximally entangled state in a <i-dimcnsional Hilbert 
space. 

In view of the recent interest in quantum information 
processing with continuous variables ||, 0, ||, ||, [To) , it is 
highly desirable to establish experimentally feasible en- 
tanglement distillation and concentration protocols for 



the continuous variables. Of particular importance are 
the protocols for Gaussian states, because these states 
can be prepared in the lab with the use of commonly 
available resources comprising lasers, passive linear op- 
tics and squeezers (parametric amplifiers). However, it 
was proved recently that it is impossible to distill Gaus- 
sian entangled states by means of Gaussian operations 
only This means that additional resources beyond 
the linear optics and homodyne detectors are required. 
The distillation protocols for Gaussian states proposed 
so far employ the photon- number measurements. The 
scheme suggested by Duan et al. relies on nondemo- 
lition measurement of the total photon number in two 
(or more) modes and represents a direct extension of 
the Schmidt projection method to infinite-dimensional 
Hilbert space. The Procrustean schemes considered by 
Opatrny et al. p"3[ ] and further analyzed by Cochrane et 
al. |Q are based on a controlled addition and subtrac- 
tion of photons. We also note that several distillation 
schemes for entangled coherent states have been proposed 

ill, n. 

In this paper, we design two entanglement- 
concentration setups for a single copy of pure two-mode 
squeezed vacuum state 



n=0 



c n \n,n) 



c n = y/l - A 2 A r> 



(2) 



where A = tanhr and r is the squeezing constant. This 
state can be generated in the process of a non-degenerate 
spontaneous parametric downconversion and provides a 
common source of the continuous- variable entanglement 
in the experiments. The Procrustean procedures that we 
are proposing preserve the structure of the state @ while 
the Schmidt coefficients c„ are transformed to new ones, 
c„ —* d n . The first scheme is based on a dispersive in- 
teraction of a two-level atom with the microwave cavity 
field and the atomic-state detection. The second scheme 
utilizes a cross Kerr interaction, coherent states, homo- 
dyne measurements, and linear optics. The underlying 
mechanism of both these schemes is that a certain aux- 
iliary system experiences a phase shift that depends on 
the number of photons in the Alice's mode of the shared 
state (|^). We convert this phase modulation into am- 
plitude modulation via interference, which allows us to 
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FIG. 1: Schematic of entanglement concentration setup in 
cavity QED: O is atomic oven, V is atomic velocity selector, 
L is laser excitation mechanism, Ri and R2 are the Ramsey 
zones driven by the microwave source, CAVITY contains the 
Alice's (Bob's) part of the entangled state and D e , D 9 are the 
field ionization detectors measuring the state of the Rydberg 
atom. 



control the amplitude of the Schmidt coefficients c n . An 
essential part of our probabilistic protocols is the mea- 
surement on the auxiliary system which tells us whether 
the concentration succeeded or failed. 

The paper is organized as follows. The first scheme is 
analyzed in Sec. II and the second scheme is discussed 
in Sec. III. Finally, Sec. IV contains the conclusions. 



II. ENTANGLEMENT CONCENTRATION IN 
CAVITY QED 
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FIG. 2: The performance of the entanglement concentration 
scheme shown in Fig. 1 for A = 1/2 and tp — n/10. (a) 
Probability of success P (solid line) and the von Neumann 
entropy 5* after the concentration (dashed line) and (b) the 
fidelity F of teleportation of coherent states are plotted as 
functions of ipo- For the input state, Si n = 0.75 and Fi n = 
0.75. 
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Our first entanglement concentration scheme is de- 
signed for the quantum state of electromagnetic field con- 
fined in a high-Q cavity and is schematically sketched in 
Fig. 1. Note, that this setup has been successfully real- 
ized experimentally and employed for the QND measure- 
ments of the cavity-field photon number and the prepara- 
tion of Schrodinger cat states [jlT], [l8| . The scheme shown 
in Fig. 1 is based on an off-resonant interaction of an (ef- 
fectively) two-level Rydberg atom with a single mode of 
a cavity sandwiched in the Ramsey interferometer. The 
atoms are emitted from an oven, their velocity is selected 
by the velocity selector, and are excited by a laser pulse 
to the upper level. Subsequently, each atom enters the 
first microwave Ramsey zone where it is prepared in a 
coherent superposition of the two (relevant) long-living 
circular Rydberg states \g) and |e), 

\<P) = ±={\g)+e^\e)). (3) 

The atom then traverses the cavity that contains the Al- 
ice's part of the shared two- mode state (||). The disper- 
sive atom-field interaction in the cavity is governed by 
the following effective Hamiltonian 

H = hna)a®\e){e\, (4) 



where a is annihilation operator of Alice's mode and k 
is effective atom- field interaction constant. The coupling 
(||) results in a phase shift Aip — a) cup of the state |e) that 
is linearly proportional to the number of photons in the 
mode A. On the other hand, the state \g) is not changed 
by the interaction. The single-photon phase shift (p = nt, 
where t is an effective interaction time, can be adjusted 
to the required value by a proper selection of the atomic 
velocity. 

After leaving the cavity, the atom passes through the 
second microwave Ramsey zone, where it undergoes a 
ir/2 Rabi rotation, 

\g)^^=(\g) + \e)), | e )^i=(|e) -|<7»- (5) 

The resulting state of the atom and the two-mode field 
reads 

00 

\^) = ^Y. c ^ l ~ elV0 ' m ' P )\9)®\n,n) 

1 00 

+ 2 E + e ltp0 ~ mtp )\e) ® \n,n). (6) 

n=0 



To complete the procedure, we measure the state of the 
atom by means of state-selective ionization detectors ex- 
hibiting almost unit detection efficiency. The entangle- 
ment concentration succeeds only if the atom is found to 
be in the ground state \g). The new Schmidt coefficients 
after this conditional transformation read 



ic n exp 



tpo - nip 



np - if 



(7) 



The irrelevant overall phase factor iexp(iip /2) can be 
dropped. Moreover, the phase factor exp(— inip/2) can 
easily be compensated by appropriate phase shift or sim- 
ply by properly redefining the quadratures of the Alice's 
mode. After these transformations, the new Schmidt co- 
efficients become real and after renormalization we get 



where 



P =2- 



-A 2 / nif _ yo 



1 1 — A 2 cos{p ) — A 2 cos{ip + cpo) 



1 - 2A 2 cosO) + A 4 



(8) 



(9) 



is the probability of success of the conditional transfor- 
mation. Clearly, two trends are competing in Eq. (^). 
The exponential decay A™ is for certain n partially com- 
pensated by the second term sm[(mp—ipo)/2] which grows 
with n up to n max = (jr + po)/p. This allows us to in- 
crease the entanglement of the shared state. 

Formally, the conditional transformation can be de- 
scribed as a diagonal filter applied to the input two-mode 
density matrix pab- Define operator 
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The output (unnormalizcd) density matrix is given by 

Pout = A ® IbpabA 1 * ® 1 B , (11) 

where Is stands for an identity operator on the Hilbcrt 
space of the Bob's mode. 

Since the conditional transformation preserves the pu- 
rity of the two-mode state, we can conveniently quantify 
the entanglement as the von Neumann entropy of the 
reduced density matrix of the Alice's mode, 
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FIG. 3: Optical implementation of the entanglement concen- 
tration scheme shown in Fig. 1. 



For the data in Fig. 2, we obtain Si n = 0.75. The fig- 
ure 2(a) clearly shows that for certain interval of phase 
shifts ipo our procedure allows us to conditionally increase 
the amount of entanglement in the pure two-mode state 
shared by Alice and Bob. 

Let us now demonstrate that the entanglement con- 
centrated in this way is useful in practical tasks. To be 
specific, we consider the teleportation of coherent states 
in the Braunstein-Kimble scheme Q where our state is 
used as the quantum channel. Making use of the transfer 
operator formalism jl9], ^0|, we can express the fidelity 
of teleportation as follows, 



F = 



1 oo oo 
m— 7i— 



m + n > 

n J 2 m+n ' 



(14) 



On inserting the Schmidt coefficients (g) into Eq. (jlj) 
and carrying out the summations we obtain analytical 
formula for the fidelity of teleportation of coherent states, 



F 



1-A 2 
AP 



1 



cos(i^o) — A cos(iy9/2 + tpo) 
1 - Acos(<ys/2) 1 - 2A 003(95/2) + A 2 

(15)" 

The fidelity F is plotted in Fig. 2(b). For fixed A and 
p we can optimize the phase po so that the teleporta- 
tion fidelity F will be maximized. For the data used in 
Fig. 2, we find that it is optimum to set po f=a — tt/10, 
which yields the fidelity F — 0.837 and the probability 
of success is P = 0.05. This should be compared with 
the fidelity F ln = 0.75 that is achieved when the original 
two-mode squeezed vacuum with A = 1/2 serves as the 
quantum channel. This improvement in fidelity is quite 
significant and clearly illustrates the practical utility of 
our procedure. 



S = -]TK|lnK| 



(12) 



n=0 



The entropy S is plotted in Fig. 2(a) as a function of ipo 
for fixed ip and A. Before concentration, Alice and Bob 
share two-mode squeezed vacuum (^|), and the entropy 
(H) reads 



S 



ln(l-A 2 ) 



1-A 2 



In A 2 . 



(13) 



III. ENTANGLEMENT CONCENTRATION FOR 
TRAVELING LIGHT FIELDS 

To implement the scheme discussed in the preceding 
Section for traveling light fields, we could replace the 
atomic Ramsey interferometer with a Mach-Zhender in- 
terferometer for a single photon and couple this auxil- 
iary photon to the Alice's mode via nonlinear medium 
with cross Kerr effect, see Fig. 3. A similar setup has 
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FIG. 4: Schematic of the entanglement concentration setup 
for traveling light fields that is based on auxiliary coherent 
states, cross Kerr interaction, eight-port homodyne detection 
(EHD), and a linear phase shift depending on the outcome of 
the measurement (PS). 



been proposed by Gerry for the generation of Schrodinger 
cat states ^J. However, this scheme has several draw- 
backs. First, the currently achievable Kerr nonlincarities 
are rather low. Secondly, we have to prepare a single 
photon. Therefore we propose an alternative scheme, see 
Fig. 4. In that setup, an auxiliary mode C is prepared 
in a (strong) coherent state |a) and then interacts with 
the Alice's mode A in the Kerr medium described by the 
Hamiltonian 



-ffKcrr = hna^ac'c. 



(16) 



After the interaction, we project the output state into 
coherent state |/3) in the eight-port homodyne detector. 

The principle of the operation of this scheme may be 
explained as follows. If there are n photons in the mode 
A, then the coherent state |a) evolves to \ae mip ), where 
ip = — Kt and t is the effective interaction time. The 
probability of projecting into |/3) is given by 



P((3\n)^-\((3\ae^}\ 2 . 

7T 



(17) 



From this formula we can see that the probability may 
grow with n if /3 belongs to certain region of the phase 
space. Without loss of generality, we may assume that 
a is real and positive and define (3 — |/3| e^q>(iipa). After 
projecting into /3, the new Schmidt coefficients can be 
expressed as 



d n oc c n ((3\ae mv ). 



(18) 



Making use of the formula for the scalar product of two 
coherent states fl23l 



we obtain 



where 



d n oc c„exp(g„ + i</>„), 



q n = | a/3 1 cos(n^ - tp ), 
4> n = | a/3 1 sin (nip - ip ). 



(19) 



(20) 



(21) 
(22) 




FIG. 5: (a) The Q-function Q([3) of the output state of the 
auxiliary mode and (b) the overlap J-(/3) are shown for A = 
1/2, q = 10, ip = 7r/100, and N = 10. The coordinates x and 
y are defined as x + iy = f3 — a. 



The old coefficients c„ are amplified or de-amplified by 
the factor exp(g„). The highest enhancement occurs for 
n = fo/ip, when q n = |a/3|. For strong auxiliary signal 
a 3> 1, the phase ipo between /3 and a will typically be 
of the order \(p Q \ w l/|a| and also |/3| w |a| will hold. 
Since the nonlinear phase shift nip will typically be very 
small for all n for which c„ substantially differs from zero, 
np <C 1, we can expand the expressions (|2l]) and ( p2[ ) in 
Taylor series and keep only terms up to linear in nip, 

q n = | a/3 1 cos(ipo) + nip\a[i\ sin(co ), 

<f>„ = -\a(3\ sin(ip ) + ntp\af3\cos(<p ). (23) 

This approximation is valid when 

\a\rup < 1 (24) 

holds. Within this approximation the exponents q n are 
linearly proportional to n and the conditional transfor- 
mation preserves the structure of the two-mode squeezed 
state: 



K| ex A" 



A = A exp(iy9|a/3| sin ip ). (25) 



Since |/3| sin cp 1, it is the product <p\a\ that de- 
termines the modulation of the input Schmidt coeffi- 
cients. A weak Kerr nonlinearity (small phase shift ip) 
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FIG. 6: The average fidelity (T) (solid line) and the probabil- 
ity of success (dashed line). All the parameters are the same 
as in Fig. 5. 



choose some reasonable figure of merit T((3) that has to 
be evaluated for each (3 and then specify a region O in the 
phase space where this figure of merit is sufficiently large. 
The entanglement concentration succeeds only if [3 € f2 
and fails otherwise. It follows that the concentration will 
yield a mixture of the states 



n, n) 



(28) 



d=0 



and the density matrix of the output state shared by 
Alice and Bob can be expressed as follows, 

Pn = 4~ [ d 2 PQ((3)mP))W)l (29) 

Here 



can be compensated by using a sufficiently strong auxil- 
iary coherent state \a). Furthermore, the strength of the 
cross-Kerr interaction can be enhanced by many orders 
of magnitude in a coherently prepared resonant atomic 
medium. A medium with electromagnetically induced 
transparency can exhibit an extremely large Kerr nonlin- 
earity g ||, ||, ^ |27| that would suffice for the practi- 
cal implementation of the present entanglement concen- 
tration scheme. 

One undesirable effect of the projection into coherent 
state is the phase modulation <p n of the Schmidt co- 
efficients. Note, however, that if the approximation j2^ ) 
holds then it follows from Eq. (|2^) that the conditional 
phase shift <f) n is linearly proportional to n and can thus 
be removed by a suitable phase shifter PS. The actual 
phase shift is proportional to the real part of (3 and we 
must use a feedforward scheme, where the operation of 
the PS (e.g., a Pockels cell) is controlled by the mea- 
surement outcome, as is schematically indicated in Fig. 
4. 

After this qualitative discussion, let us now provide a 
rigorous mathematical description of the entanglement 
concentration scheme shown in Fig. 4. The probabil- 
ity that (3 will be measured in the EHD is the sum of 
all the conditional probabilities P((3\n) multiplied by the 



probabilities 



c n \ that there are n photons in the mode 



Q(P) 



i-x 2 



^A 2 "exp(-|ae 



inip 



(26) 



n=0 



The normalized Schmidt coefficients corresponding to the 
measurement outcome [3 read 



dn(J3) 



VI - A 2 A™ exp [a(3*e mv - inip\a(3\ cos ip ] 
V^Qp)exp(H 2 /2-H/?| 2 /2) 



(27) 

We need to establish a criterion according to which we 
will accept or reject the state in dependence on the mea- 
surement outcome (3. The most natural approach is to 



Pn = 



d 2 {3Q{(3) 



(30) 



denotes the probability of success of the concentration, 
i.e., the probability that the measurement outcome f3 will 
belong to il. 

Different figures of merit may be suitable depend- 
ing on the intended usage of the shared quantum state. 
For instance, if that state shall be used to teleport coher- 
ent states, then it will be natural to employ the fidelity 
@ as the fi gure of merit. Here we make use of an even 
simpler quantity, namely the fidelity 



^(/?) = i<<«(/3)>r 



(31) 



between the conditionally prepared state and a maxi- 
mally entangled state in the Hilbert space of the first 
N + 1 Fock states, 



1 



n=0 



n, n) 



(32) 



The fidelity depends on the sum of the first A^+l Schmidt 
coefficients and is explicitly given by the following for- 
mula, 



1 



N + 1 



N 



n=0 



(33) 



In the rest of this section we present the results of 
the numerical calculations for a — 10, tp = 7r/100, and 
A = 1/2. The Q-function (|2^) of the output auxiliary 
mode is plotted in Fig. 5(a). Since the assumed nonlinear 
single-photon phase shift 7r/100 is relatively small, the 
Q-function is practically identical to the input Gaussian 
Q in (/3) = exp(-|a-/3| 2 )/7r. The function is plotted 
in Fig. 5(b) for N = 10. We can see that there are 
regions in the phase space where T is higher than the 
fidelity corresponding to the input two-mode squeezed 
vacuum, Tq = 0.273. As described above, we define 17 as 
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FIG. 7: The fidelity of the teleportation of coherent states 
when the state pn after concentration is used as the quantum 
channel. All parameters are the same as in Fig. 5. 

the region of the phase space where F((3) > Fo + AF. 
The dependence of the average fidelity 

(f) = TT f d 2 (3Q(P)F((3) (34) 

on AF is plotted in Fig. 6. Also the probability of suc- 
cessful entanglement concentration ( |30| ) is shown there. 
As the gap AT becomes larger, the average fidelity in- 
creases while the probability decreases. 

Finally, we show that the entanglement concentrated 
in this way is suitable for the teleportation. We can cal- 
culate the average teleportation fidelity similarly as {F) . 
We evaluate the fidelity of teleportation F{j3) for each 
(3 E £1 by inserting the relevant Schmidt coefficients d(j3) 
given by Eq. (p7j) into Eq. and then we average F(J3) 
over with the properly normalized probability density 



Q{P)/Pq. The results are shown in Fig. 7. We can see 
that the teleportation fidelity monotonically grows with 
AF and for all AF > we have F > 0.75. This ex- 
ample confirms that our procedure indeed extracts more 
useful entanglement from the input two-mode squeezed 
vacuum. 



IV. CONCLUSIONS 

In this paper, we have designed two schemes for prob- 
abilistic concentration of continuous-variable entangle- 
ment. The Procrustean protocols that we are proposing 
have the important property that they can be applied 
several times to a single copy of the shared two-mode 
entangled state. Thus we could, in principle, extract a 
state with very high entanglement, at the expense of a 
low probability of success. When repeating the concen- 
tration procedure, one could optimize the relevant pa- 
rameters such as the phase shifts (fo and ip in order to 
achieve the optimum performance of the schemes. In 
view of the recent advances in cavity-QED experiments 
and the preparation of media with extremely high Kerr 
nonlinearity, we may hope that the schemes proposed in 
the present paper will become experimentally feasible in 
a near future. 
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